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Abstract 

The Clifford algebra of a n-dimensional Euclidean vector space provides a 
general language comprising vectors, complex numbers, quaternions, Grass- 
man algebra, Pauli and Dirac matrices. In this work, a package for Clifford 
algebra calculations for the computer algebra program Mathematicd^ is in- 
troduced through a presentation of the main ideas of Clifford algebras and 
illustrative examples. This package can be a useful computational tool since 
allows the manipulation of all these mathematical objects. It also includes 
the possibility of visualize elements of a Clifford algebra in the 3 -dimensional 
space. 
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1. Introduction 



The importance of Clifford algebra was recognized for the first time in 
quantum field theory but there has recently been a tendency to exploit their 



power in many others fields. These fields include projective geometry 13], 
electrodynamics |l4] • analysis on manifolds and differential geometry 12], 
crystallography [2j, |3] and others. A recent account on applications of Clif- 
ford algebra in fields such as robotics, computer vision, computer graphics, 
engineering, neural and quantum computing, etc., can be found in [5j and 
3- 



Libaries, packages and specialized programs for doing Clifford algebra 
already exist; CLIFFORD /Bigebra is a Maple package which includes ad- 
ditional specialized packages such as SchurFkt (for the Hopf Algebra of 
Symmetric Functions) and GfG - Groebner for Grassmann (for Computing 
Groebner Bases for Ideals in Grassmann Algebra) [1]; I "Clifford Algebra is an 
add-on application for the Tensorial Mathematica package that implements 
some Clifford algebra operations [8J Clifford is a Mathematica package for 
doing Clifford algebra, oriented to differential geometry [l8| and CLICAL is 



a stand-alone calculator-type computer program for MS-DOS [17 . 

While the first two packages requires more specialized knowledge on Clif- 
ford algebras, CLICAL is easy to use and can be used by non mathematicians. 
It has the advantage of being a stand-alone program (does not require a pro- 
prietary program to run) but does not allow symbolic computations and, 
as far as we know, is currently not under maintaining. The Mathematica 
package Clifford described in |18j (which is not easily available) performs 
symbolic calculations and multivectors notation is quite similar to the used 
in our package. The main difference lies in the implementation; while the 



package by Parra and Rosello [18] is based on matrix representations of the 
Clifford algebra, cliff ord.m is based on isomorphisms between blades and 
n-tuples. As we shall see, this approach enables a more efficient and elegant 
computational scheme, providing also further flexibility concerning the de- 
pendence on the dimension of the space under work. In most of the package 
functions, for instance, the dimension of the space where a multivector is 
embedded is calculated on the fly, so it is not necessary to specify since the 
very beginning the dimension of the space under work. 

The Mathematica package herein described has the main purpose to im- 
plement general operations of a Clifford algebra on the language of the com- 
puter algebra program Mathematica. In order to introduce the package and 
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its main algorithms, we first present a gently introduction to the Clifford al- 
gebra of M n , with examples thought to show the generality of this algebra in 
the sense of providing a general language comprising vectors, complex num- 
bers, quaternions, Pauli and Dirac matrices. These examples are worked out 
with cliff ord.m with the aim of describing the use of the main predefined 
functions. For more general introductions to Clifford algebras, readers are 



referred to books, such as Refs. [19( and [9|] 

We must point out that an early version of cliff ord.m, intended for ap- 
plications in crystallography, was given in Ref. Q. The package presented 
here, however, is more general in purpose; the code has been rewritten and it 
has been enriched with functions to draw multivectors in M 3 . The program 
has been tested in several ways. For instance, calculations and graphical 
examples in the work by Zhang et al. ji^] were repeated and the resulting 
notebook is available for download (see below). The package has been also 
used in works by our group 0, 0, 0] and also in Ref. (24f . The package 
cliff ord.m, a user guide (UserGuide .nb), a Palette with the most com- 
mon predefined functions (Clif f ord.nb) as well as the notebook with the 
calculations by Zhang et al. |25( (ZhangEtAl .nb) are available for download 
at: 

\protect\vrule widthOpt\protect\href -fhttp : //www. fata. unam.mx/ar agon/ sof tware/Mhtt 



2. The Clifford algebra of R" 

Let us consider the vector space MJ 1 with the inner product (, ) and an 
orthonormal basis {ei, e2, e n }. One of the most simple and direct ways to 
define a Clifford Algebra is by making use of their generators [g]. 

Definition 2.1. The real associative and distributive algebra generated by 
the Euclidean space M. n with the product rules 

e\ = 1, * = l,2,...,n, (1) 
a^Gj -)- e^e^ 0, % 7^ j, 

where {ei, e 2 , . . . , e n } is the canonical basis ofW 1 , is called universal Clifford 
algebra of the space M. n and is denoted by M. n . 

With this definition, we can easily establishes the relationship between 
the canonical scalar product in MJ 1 and the so-called geometric product of 
the algebra R n . 
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Proposition 2.2. I/x, y e R n i/ten x > and xy + yx £ R. Even more 

x 2 = (x, x) 

xy + yx 

— 2"" = (x ' y) 
where (■, ■) is the canonical inner product in R n . 

Proof. It is enough to write x and y in terms of the canonical basis 
{ei, e 2 , . . . , e n } and use Equation (pQ). 



In a Clifford algebra the inverse of a nonzero vector can be defined [6j. 
In particular, the inverse of any vector x 6 R™, x ^ 0, is given by 



x x , n \ 

X - := = (2) 



1 _ 



,, X ,, 

The Clifford algebra R n is itself a vector space of dimension Ylp=o («) 
2 n , with basis 

{1, ei, e 2 , . . . , e n , eie 2 , eie 3 , . . . , eie n , . . . , • • • e n } , 
such that an element A in R„ is written as 



A = a + anei H h a u e n + a 2 ieie 2 H a 2i e 1 e n H h a di eie 2 • • • e n , (3) 

where d = 2 n — 1 and z = (™) for the real numbers a V {. Consequently, the 
vector space R n can be decomposed, as direct sum, in n + 1 subspaces: 

R n = A°R n © A x R n © • • • © A n R n . (4) 

where A corresponding to the exterior product 0] and each subspace is of 
dimension (™). 

The elements A (Eqn. [3D of the Clifford algebra R n are called multivec- 
tors, and those of A p R n , p-vectors. In particular, 0-vectors are real numbers 
and dim(A°R n ) = 1. A 1 R ra has the basis {ei, e 2 , . . . , e n }, so 1-vectors are sim- 
ply vectors and dim(A 1 R n ) = n. A 2 R n has the basis {eie 2 , eie 3 , . . . , eie n } 
and their elements (2-vectors) are also called bivectors. Finally, A n R n has 
as basis {eie 2 ■ • • e n } and since dim(A n R n ) = 1, the n-vectors of R n are also 
referred as pseudoscalars. 
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In what follows, arbitrary multivectors will be denoted by non bold upper- 
case characters without ornamentation such as A. p-vectors will be denoted 
by A p , with the exception of vectors (1-vectors), that will be denoted by bold 
lowercase characters such as a. 

Bearing in mind the decomposition of the vector space M. n as the direct 
sum of the subspaces A P R™, < p < n, given in (Jl}, any multivector A can 
be written as 

A=(A) + (A) 1 + --- + (A) n , (5) 

where (A) p , the p-vector part of A, is the projection of A e M. n into A p IR n . 
( ) is called the grade operator. 

2.1. Inner and outer products 

Given the decomposition (j3J), an important property of a Clifford algebra 
is the existence of products that allows us to move from one subspace of 
M. n to another. Let us first consider the product of two 1-vectors. For all 
u, v G A 1 !^, their product uv can be written as 

uv = - (uv + vu) + - (uv — vu) . 

Now define the "inner" and "outer" products as follows 

u ■ v = - (uv + vu) = (a, b), 

1 , 

uAv = - (uv — vu) . 

The inner product is symmetric and notice that vectors u and v are orthog- 
onal if an only if uv = — vu. The outer product u A v is antisymmetric (and 
associative) and vanishes whenever the two vectors are collinear, that is, u 
and v are collinear (or linearly dependent) if an only if uv = vu. Thus, 
the product uv provides information about the relative directions of the 
vectors. Anticommutativity means orthogonality and commutativity means 
collinearity. Notice that for the basis vectors of M n , we have 

6j6j Gj A Gj, % ^ J, 

From the following equality 

uv = u ■ v + u A v, 

= (UV) + (UV> 2 , 
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we can extend the notions of inner and outer product to the case of p- and 
g-vectors in the following way. For a p-vector A p e A p M. n and a g-vector 
B q e A 9 IR n , the inner product A p ■ B q is defined by 

f (A v B a ), | if p,(/>0, 
p 9 ~ \ if p = or g = 0. 1 ' 

Analogously, the outer product A p A B q is defined by 

A p A B q = (A p B q ) p+q . (7) 

Since arbitrary multivectors can be decomposed as in (jSJ), inner and outer 
product can be extended by linearity to R n . Then, given A,B e M n , we have 

n 

A.S=53<A> Jfc .(S> 1 , (8) 

fc,Z=l 

?t 

AAfi= J2( A )^( B )h (9) 

fc,i=i 

S.S. Geometric interpretation 

Bivectors have an interesting geometric interpretation. Just as a vector 
describes an oriented line segment, with the direction of the vector repre- 
sented the oriented line and the magnitude of the vector is equal to the 
length of the segment; so a bivector a A b describes an oriented plane seg- 
ment, with the direction of the bivector represented the oriented plane and 
the magnitude of the bivector measuring the area of the plane segment (Fig- 
ure [TJ. The same interpretation is extended to high-order terms: a A b A c 
represents an oriented volume. In general, multivectors contain information 
about orientation of subspaces. 

2.3. The reverse of a multivector 

The magnitude or modulus of a multivector A is defined by the equation 

/ ~ \ 1/2 

\A\ = (AA) Q , (10) 

where ~ denotes the operation reverse defined as 

(eie 2 - ••e p )~ = e p - ■ ■e 2 e 1 . 



The operation reverse is distributive [12[ so that the reverse of an arbitrary 
multivector A can be easily calculated. 

If the inverse of a multivector A exists, it is denoted by A^ 1 or 1/A, and 
is defined by the equation AA^ 1 = 1. 
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2.4- General metrics 

In many applications one considers real vector spaces W 1 with metrics 
that are not positive definite with the bilinear form ( , ), such that 

/ \ — 2 _i_ 2 _i_ _i_ 2 _ 2 _ _ 2 

where n = p + q. In this case, the vector space is denoted as W' q , giving rise 
to the Clifford algebra K Pi(J |6| : 

Definition 2.3. The real associative and distributive algebra generated by 
the space M. p,q (p + q = n) with the product rules 

e- = 1, i = l,2,...,p, 

e 2 = -1, i = l,2,...,q (11) 
+ e^e, = 0, i ^ j,i, j = 1,2, ... ,n 

where {ei, e 2 , . . . , e n } zs £/ie canonical basis ofW' q , is called universal Clifford 
algebra of the space W' q and is denoted by R Pi9 . 

As in Proposition 12.21 the relationship between the canonical scalar prod- 
uct in M. p ' q and the geometric product of the Clifford algebra can be 
established. 

Proposition 2.4. /fx, y G M. p,q then x 2 and xy + yx G K. i£ven more 

x 2 = (x, x) 
x y + ^ x / \ f-to\ 

— 2 — = ( x >y)> ( 12 ) 

where (•, •) is the canonical inner product in M. p,q . 

Proof. If x = Y^l=i x i e i anc ^ y = J2i=i Vi e i then, by ffTTj) we have 

P p+g 

x 2 = Y, x t- E x * 2 e M 

i=2 i=p+l 

p p+g 
xy + yx = Sij/i - ar^ G R. 

i=2 i=p+l 

Now, an element a G M p ' 9 is said to be invertible if a 2 ^ 0. The element 

a 



-l _ 

o2 



is called the inverse of a. 
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3. Clifford algebra calculations with Mathematica 

According to Eqns. §}, (JSJ, ®, ©, © and all that we should 

need in order to manipulate multivectors in a computer algebra program 
such as Mathematica, would be to define the two fundamentals operations 
(algorithms): geometric product and grade operator. 

3.1. Fundamentals algorithms 

The first fundamental algorithm for the computation of the geometric 
product between multivectors can be devised by noticing that a general mul- 
tivector A in M. PtQ , such as (J3j), is formed by a linear combination of terms of 
the form 

eJ*eS*...c£-. (13) 

where = 1,0, (i = 1, n), that we call blades. The geometric product of 
two of these blades (which is also a blade) is: 



(ePeS* • ■ ■ e-r ) (e?e r 2 2 ■ • • <») = (-l^ep^eS" 4 " ■ ■ ■ C^", (14) 
where the sum m; + r» is evaluated modulus two, and 



S = > TiVTLy 



l<i<j<n 

If rrii + ?"j = 2 then, in order to have the right hand side in the form ( {TBI 
when considering the signature of the bilinear form, e™ i+Ti will be replaced 
with (ej, e,)e°, and in this case we have: 

mi+n m,+r, „m n +r n _ / „ \ c> m i+ r 'i o „m. n +r n 

e l ' ' ' e i ' ' ' e n — \ e i> " ' e ! " ' e n 

Equation (fT3|) enables us to establish an isomorphism between blades and 
n-tuples 

e^e™*---e^^( mi ,m 2 ,...,m n ), 

that can be more easily manipulated from a computational point of view, 
represented by lists of S and l s . 

As an example, let A 2 = 7eie 2 and B 2 = 636-4 be two blades in M4. 
They can be represented as e^e^e^e^, (1,1,0,0), and e^e^e^e^, (0,0,1,1), 
respectively. The geometric product A 2 B 2 is given by: 

A 2 B 2 = (-I) 2 7ej +1 e0 +1 ei +0 e0 +1 = 7e 2 e 3 e 4 , 
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or (0,1,1,1). 

Once we are able to calculate the geometric product between two arbitrary 
blades A p and B q in K n , the product between two multivectors A and B in 
R n can be evaluated by using the decomposition (jSJ). Indeed, if 

?! 

i=0 
n 

3=0 

the geometric product AB is given by 

n n / n 

AB = j2J2(- i r^ IK 

8=0 j=0 \fc=l 

where s — ^2 rj p rrii q and m iq + r Jp = mod2. 

l<p<g<n 

Finally, according to (jSJ), ©, and ffTUl) . to calculate the inner and outer 
product between general multivectors and the magnitude of a multivector, 
the grade operator must be implemented (also the reverse operation but it 
turns to be trivial). Given a multivector A £ Mp g , to calculate, for instance 
(v4) r , we proceed as follows 

1. Decompose A into a sum of blades as in ((Sj). 

2. According (fl3l) . the grade of the i-th blade of the above decomposition 
is simply run + m i2 + • • • + m in . 

3. Keep the terms of the decomposition with grade equal to r. If no such 
terms exist then (A) r =0. 

3.2. Implementation 

We denote, in Mathematica code, the j-th basis vector as e [j] , a blade 
such as e!e 3 e4, ( e^e^ege^ using the nomenclature of Eqn. [TBI is written as 

e[l]e[3]e[4], (15) 

and can be internally represented simply by (1, 0, 1, 1). Care must be taken in 
preserving the canonical order of the expression, since for instance e [1] e [3] 
is a geometric product of two vectors and e [1] e [3] ^ e [3] e [1] . 



i=0 \fc=l 



k 



Eft IK * 

j=0 \fc=l 
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Let us consider the Clifford algebra M P)9 . If dim = p + q is the dimen- 
sion of the vector space, the following Mathematica code implements the 
transformation of a blade, such as (fT5j) . onto a n-tuple: 

ntuple[x_, dim_]:= ReplacePart [Table [0, {dim}] , 1, List @@ x /. e[k_] -> -Tk}] 

The signature of the bilinear form can be set by using $SetSignature = p 
and Equations fllip can be coded as: 

bilinearf orm[e [i_] ,e [i_] ] := If [i <= $SetSignature , 1, -1] 

Except for some functions, it is not necessary to define the dimension of 
the vector space since it can be calculated directly The maximum dimension 
of the space where a blade is embedded can be extracted from the list: 

dimensions [x_] : = List @@ x /. e [k_?Positive] -> k 

To calculate the dimension of the space where a a general multivector 
is embedded, the previous function must be generalized by including the 
distributivity of addition: 

dimensions [x_Plus] : = List @@ Distribute [f [x] ] /. f -> dimensions 

dimensions [a_] : = {0} /; FreeQ [a,e [_?Positive] ] 

dimensions [a_ x_] := dimensions [x] /; FreeQ [a, e [_?Positive] ] 

3.2.1. The geometric product 

From (114j) . the geometric product between two blades with a bilinear form 
of signature p can therefore be evaluated with help of the function geoprod: 

geoprod[x_,y_] := Module [{q=l, s, r=0, 

pl=ntuple[x, Max [dimensions [x] , dimensions [y] ]] , 
p2=ntuple[y, Max [dimensions [x] , dimensions [y] ]]} , 
s=Sum[p2[[m]]*pl[[n]] , {m .Length [pi] -1}, 

-[n,m+l, Length [p2]}] ; 

r=Mod[pl+p2,2] ; 

Do [If [rl [ [i] ] == 2 , q*=bilinearf orm [e [i] , e [i] ] ] ; 
If [r[[i]]==l, q*=e[i], {i, Length [r] }] 
(-l)-s*q ] 
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This function evaluates the geometric product of two blades in One 
step further consists in to calculate the geometric product of two arbitrary 
multivectors such as ([3]). This can be achieved from geoprod by providing the 
transformation rules which contains the properties of the geometric product 
under multiplication of blades by real numbers and addition of blades. Here 
is the behavior under scalar multiplication: 

geoprod [a_, y_] := ay /; FreeQ [a, e [_?Positive] ] 

geoprod [x_, a_] := ax /; FreeQ [a, e [_?Positive] ] 

geoprod[a_ x_,y_] := a geoprod[x,y] /; FreeQ [a, e [_?Positive] ] 

geoprod[x_,a_ y_] := a geoprod[x,y] /; FreeQ [a, e [_?Positive] ] 

and the distributivity of addition: 

geoprod [x_, y_Plus] := Distribute [f [x,y,p] , Plus] /. f -> geoprod 
geoprod [x_Plus, y_] := Distribute [f [x,y,p] , Plus] /. f -> geoprod 

Therefore, the function to calculate the geometric product of arbitrary 
multivectors is defined as 

GeometricProduct [_] : = $Failed 

GeometricProduct [ml_,m2_,m3 ] := f[ GeometricProduct [ml ,m2] ,m3] 

/. f -> GeometricProduct 
GeometricProduct [ml_,m2_] := geoprod [Expand [ml] , Expand[m2]] 

3.2.2. The grade operator 

To complete the basic operations of a Clifford algebra, it remains to im- 
plement the grade operator. The following auxiliary function calculates the 
grade of a blade in E p>g 

gradblade [a_] : = /; FreeQ [a,e [_?Positive] ] 
gradblade [x_] : = Plus @@ ntuple[x, Max [dimensions [x] ] ] 
gradblade [a_ x_] := gradblade [x] /; FreeQ [a, e [_?Positive] ] 

Now, Grade [x , n] should extract the term of grade n from the multivector 
X. Firstly, we consider the case when the multivector X is a blade of grade, 
say r: (X) n = if r ^ n and (X) n = X if r = n. The code reads: 

Grade [x_, n_?NumberQ] := If [gradblade [x]==n,x,0] ] 

For a general multivector X, we have: 
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Grade [x.Plus , n_?NumberQ] := Distribute [f [x,n] , Plus] /. f -> Grade 

Functions GeometricProduct and Grade enable us to construct all the 
operations which can be defined in a Clifford algebra, such as outer product 
(OuterProduct [v, w, . .]), inner product (innerProduct [v, w, . .]), magni- 
tude (Magnitude [v] ), reverse (Turn [v] ), inverse (Mult ivector Inverse [v] ), 
dual (Dual [v,dim] ), and many others, all included in the package Clif f ord.m. 
This package works with general multivectors of the form ([3]), but particu- 
lar cases can help to envisage the power of multivector calculus, as is made 
explicit in what follows. 



4. Vectors 

Let us consider the Clifford algebra M Pi o- n-dimensional vectors are 1- 
vectors and lie in the subspace A 1 M P ' . A vector a is therefore a = (A)x, 
where A is a general multivector. 

The inner product defined in ([6]) becomes now the standard "dot prod- 
uct" between vectors. The "cross product" a x b of vector calculus is defined 
in M. 3 and is related to the outer product of vectors a and b as follows. 
The vector a x b is perpendicular to a A b and with the same magnitude: 



|aAb| = |axb|. The explicit algebraic relation between them is [11 



a x b = (-eie 2 e 3 )(a A b), (16) 

where e^ea is the pseudoscalar of IR^o- We may actually take this as a 
definition of the cross product. 

From (1161) it is easy to define the function crossprod [v,w] that gives the 
cross product between two three-dimensional vectors v and w: 

crossprod [v_ , w_] : = GeometricProduct [-e [1] e [2] e [3] , OuterProduct [v , w] ] 

The associativity of the geometric product allows algebraic manipulations 
typical of real numbers that are no possible in the Gibbs' vector algebra since 
the cross and dot products are not generally associative. For example 

a x (b x c) ^ (a x b) x a. 

Even more, many products such as a • (b • c) are not even defined. With 
Clifford algebra all products are not only well defined but associative mak- 
ing simpler many algebraic manipulations and allowing to define derivatives 
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and integrals just as they are defined for real functions of real variables, pro- 
vided that we are careful to maintain the order of the factors since geometric 
product is not commutative. 

One specific example that shows the simplicity of some expressions if 
geometric product is used, is concerning rotations. Consider a vector v in 
R n which is rotated by an angle 9 in the oriented plane characterized by the 
bivector a A b. After the rotation, the vector v is transformed into v', given 
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v = UvU, (17) 

where 

U = cos(fl/2) + sin(0/2). 
|aAb| 

The direction of the rotation (clockwise or counterclockwise) is specified by 
the orientation of the bivector a A b. Eqn. [17] gives the rotated vector v' 
with no matter the dimension of the space in which it is embedded. No 
corresponding simple expression exists in vector algebra. 

As a quite simple example, which can be easily visualized, consider the 
vector v = (1,1,1), to be rotated 90° maintaining invariant the plane xy. 
To characterize the plane xy, we can use a A b where a = (1,0,0) and 
b = (0, 1, 0), in which case we get a rotation counterclockwise and is easy to 
see that v' = (—1, 1, 1), but if we use b A a then the rotation is clockwise 
and v' = (1, —1, 1). Here is this example solved with Clifford: 



In[l] 
In [2] 
In [3] 



= « Clifford. m 

= v = e[l]+e[2]+e[3] ; 

= plane = OuterProduct [e [1] ,e [2] ] /Magnitude [OuterProduct [e [1] , e [2] ] ] 



Out [3]= e[l]e[2] 

The operator U is now defined (U is the reverse of U): 
In [4]:= u = Cos [Pi/4] + plane Sin [Pi/4] ; 

In[5] := vprime = GeometricProduct [Turn[u] ,GeometricProduct [v,u] ] 
Out [5]= -e[l]+e[2]+e[3] 

So we get a counterclockwise rotation where the vector v = (1, 1, 1) be- 
comes v' = (—1, 1, 1). The implementation in Mathematica of this rotation 
is the function Rotation[v,a,b] . 
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5. Drawing multivectors in the 3-dimensional space. 

The elements of the Clifford algebra l^o can be visualized in the 3- 
dimensional space, providing geometrical insights. The package include a 
function to draw multivectors belonging to l^o, called GADraw, that will be 
described in what follows. 

GADraw includes embedded functions to draw vectors, bivectors and the 
pseudoscalar of l^o- For example let us draw the multivector e 3 + eie 2 + 
e 2 e 3 + eie 2 e 3 : 

In[l] := « Clifford. m 

In [2] : = A = e [3] +e [1] e [2] +e [2] e [3] +e [1] e [2] e [3] ; 

In [3] : = GADraw [A] ; 
Out [3] = 



Scalar = 




1 



Figure 1: A vector, a bivector (plane) and the pseudoscalar (cube) drawn with the aid of 
the function GADraw. 

The result is shown in Fig. [IJ The bivector is represented by an area 
and the pseudoscalar as a scalable cube. In the particular case of the vector 
(0, 0, 1), the arrow in its tip was generated by the following code((23|): 
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mat[l] = Sin [t]*(e [1]/14) + Cos [t] * (e [2] /14) , 

mat [2] = Sin[t + 0.25] *(e [1] /14) + Cos [t + 0.25] *(e [2] /14) , 

mat [3] = e[3]/5, 

This arrow is then translated and rotated (in this case to the tip of the vector 
(0, 0, 1)) with the aid of the function Rotation [mat , w,p] , where w = e% and 
p is the vector that points the site where the tip of the arrow must be located. 
If sc = Sqrt[p[[l]] ~2 + p[[2]]~2 + p[[3]]~2]/2 is a scale factor then 
all this procedure can be encoded as: 

If [OuterProduct [ToBasis [p] , e[3]] ===0, 

cone = Table [Array [ToVector [mat [#] ,3] &,3] +p-ToVector [mat [3] ,3] , 

{t,0.25,2*Pi,0.25}] , 
elms = Array [sc*ToVector [Grade [Rotation [mat [#] ,e[3] ,ToBasis[p]] ,1] ,3]&,3] ; 
cone = Table [elms+p-res [ [3] ] , {t , . 25 , 2*Pi , . 25}] ] ; 

arrow = Graphics3D [{FaceForm [color] , EdgeForm[] , Polygon /@ cone}] 

Where ToBasis and ToVector are functions, defined on the same package 
Clifford; the former changes from the coordinates of a vector to the basis 
notation (fl5|) and the last one does the opposite. 

It is worth to mention that the above coded was included with the aim 
to show that symmetry operations, rotations in this case, can be carried 
out by using the elements of the Clifford algebra without requiring matrices. 
This algebra provides a consistent computational framework with significant 



applications in computer graphics, vision and robotics 10 



6. Complex numbers 

Let us consider the Clifford algebra of the most simple space that has a 
geometrical structure: the plane R 2,0 . Taking the canonical basis {ei,e2}, a 
basis for the Clifford algebra R2,o is {1, ei, e 2 , eie 2 } and a general multivector 
A G M2,o has the form 

A = k + kiei + A; 2 e 2 + k^eie 2 - 

We can decompose M 2j o as M.2,0 = I&2 0©^2 0> sucn ^ na ^ ^2~o contains even 
grade elements and M 2 contains odd grade elements. Therefore, A can be 
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expressed as the sum of two multivectors: A = A + + A , where A + G 
and A" e R^o- That is 

A = A + + A-, 
A + = k + fc 3 eie 2 , 
A~ = fciei + fc 2 e 2 . 

We focus our attention into ffi^o! it i s itself an algebra so that it is called 
the even subalgebra of R 2j o- By taking 

i = eie 2 , (18) 

an element of can be written as 

z = k + k 3 i. 

Since i 2 = (eje 2 )(eie 2 ) = — (eie 2 e 2 ei) = —1, and i is itself a generator 
of rotations [ll|, we see that is equivalent to the algebra of complex 
numbers. 

The algebraic operations of complex numbers can therefore be worked out 
with Clifford. In order to get a more standard notation we firstly define 
the function 

In [6]:= transform [x] := x /. i -> e[l]e[2] 

which makes the identification (1181) . Here we have two complex numbers: 

In [7] := w = a + b i; 
In [8] := z = c + d i; 

The product (a + bi) (c + di) becomes 

In[9] := GeometricProduct [transform [w] , transform [z] ] /. e[l]e[2] -> i 
Out [9]=ac-bc + bci + adi 

The equivalence between some built-in basic operations of complex num- 
bers in Mathematica and those of a Clifford algebra defined the package is 
shown in the following table: 
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Built-in Objects Cliff ord.m Objects 



Conjugate [z] 
Abs[z] 



Re[z] 
Im[z] 



Grade[z,0] 

GeometricProduct [Grade[z,2],-e[l]e[2]] 

Turnfz] 

Magnitude [z] 



The inverse of the complex number w can be calculated with Multivectorlnverse: 
In [10] := Multivectorlnverse [transform [w]] /. e[l]e[2] -> i 

a - bi 



7. Quaternions 

Taking further the same procedure, we now develop the Clifford algebra 
of M 3 ' , equipped with the canonical basis {ei, 62,63}. The Clifford algebra 
IR 3j o is an eighth-dimensional vector space with the basis 

{1, ei, e 2 , e 3 , eie 2 , e 2 e 3 , eie 3 , eie 2 e 3 }. 

A general multivector A e K.3,0 is written as 

A = k + fciei + A; 2 e 2 + k 3 e 3 + fc 4 eie 2 + k 5 e 2 e 3 + /c 6 eie 3 + k 7 e x e 2 e 3 , 

which can be also expressed as A = A + + A~ , where 

A + = k + /c 4 eie 2 + k 5 e 2 e 3 + A; 6 eie 3 
A~ = fciei + A; 2 e 2 + k 3 e 3 + /c 7 eie 2 e 3 . 

The even-grade elements A + form the subalgebra M^ of R 3;0 , equivalent 
to the algebra of quaternions. This can be seen by making the identifications 

i = -e 2 e 3 , 

j = e!e 3 (19) 
k = -eie 2 , 



Out [10] = 



2 

a + b 



2 
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leading to the famous equations 

i 2 = f = k 2 = -l, (20) 
ijk = -1. 

With the identifications given in ffl9l . an element of BU!~ can be written 
now as 

Q = Qo + <?ii + Q2i + <?3k, 
which, in view of the properties (1201) . is a quaternion. 

The algebra of quaternions is therefore comprised in the same package. 
The basic operations of this algebra are carried out by the function al- 
ready defined, such as GeometricProduct, Multivectorlnverse, Magnitude 
and Turn. To simplify the operations, we have incorporated the defini- 
tions (Tl9l) and redefined some functions to work only with quaternions and 
complex numbers. These new functions begin with the word Quaternion, 
namely, QuaternionProduct, Quaternionlnverse, QuaternionMagnitude 
and QuaternionTurn. So, for instance the inverse of the quaternion q = 
a + 3i + 6j — 10k is 

In [11] := q = a + 3 i + 6 j - 10 k; 
In [12] := Quaternionlnverse [q] 



a - 3 i - 6 j + 10 k 



Out [12] = 



145 + a 



8. Grassmann algebra 

The outer product (OuterProduct) defined in (J7]) is associative and the 
identity 1, in a Clifford algebra R^o, is also the identity for the outer product. 
Consequently, the vector space W'° with the outer product already defined 
is an algebra of dimension 2 ra , which is called the Grassmann algebra of M p '°. 
Notice that it does not depend on the inner product of the vector space, but 
just on the alternation of the outer product. 

Grassmann algebra has relevance in modern theoretical physics [16] and 
has the right structure for the theory of determinants. That the structure of 
this algebra is contained in Clifford algebra allows to reformulate Grassmann 



calculus [16L and to give an extensive treatment of determinants [12J, both 



in terms of Clifford algebra. 
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9. The hyperbolic plane 



Perhaps the simplest example of a problem involving non positive definite 
metrics is the Minkowsky model of the hyperbolic plane. Here we shall 
develop some basic ideas and calculations just to give a flavor of this kind 
of application of Clifford algebra and the use of the package Clifford for 
non positive definite metrics. In the next section, one extra dimension is 
introduced and the resulting algebras have great relevance in relativity theory 
and quantum mechanics. 

One can always visualize a surface of constant positive Gaussian curvature 



1/R 2 (see for instance 21] for concepts related to differential geometry) as a 
sphere, with radius R, embedded in a three dimensional Euclidean space. The 
surface is described by the equation x\ + x 2 + x§ = R 2 . A surface of constant 
negative curvature, however, can not be embedded in a Euclidean space, so 
alternative possibilities must be developed to visualize such surfaces. The 
simplest surface of constant negative curvature is often called the hyperbolic 
plane, the Bolyai-Lobachevsky plane, or the pseudo sphere. This surface can be 
globally embedded in a space equipped with the Minkowsky metrics instead 
of the Euclidean one. A three-dimensional Minkowsky space can be identified 
by the fact that if (x\, x 2 , £3) are the coordinates of a vector x in this space, 
then the distance to the origin is |x| 2 = x\ + x\ — x\. 
The equation 

xj + x 2 2 -xj = -R 2 , (21) 

defines a hyperboloid of two sheets intersecting the X3 axis at the points ±1. 
Either sheet (upper or lower) models an infinite surface without a boundary 
(the Minkowsky metric becomes positive definite upon it) that, as we shall 
see, has constant Gaussian curvature —1/R 2 . 

We can easily convince ourselves that M.2,i is indeed the three-dimensional 
Minkowsky space (the signature of the bilinear form is 2). We shall proceed 
to calculate the Gaussian curvature of the hyperboloid (I2T]) with the stan- 
dard formulas of differential geometry but with the metrics of IR 2 ,i • hi three 
dimensions ,_the Gaussian curvature of a surface f(xi,x 2 ,x 3 ) = can be 
written as 
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fc = [ n - V 2 n+(V-n) 2 ] , 

where n = Vf(x%, x 2 , x 3 ) / | V/(xi, x 2 , x 3 )\ is the normal to the surface. In 
a three-dimensional space with a non positive definite metric, the gradient 



19 



of a scalar function and the divergence and Laplacian of a vector function 
f = (A, A, A) are defined as 

_ , , x dcj) dcj) . . dcj) 

V0 = (e 1 .e 1 )^e 1 + (e 2 .e 2 )— e 2 + (e 3 -e 3 )— e 3 , 

V7 C I \2^A , / \2^A . / \2^A 

V-f = (ei ■ eO ^ + (e 2 • e 2 ) — + (e 3 ■ e 3 ) — , 



V 2 f = 



<9 2 B 2 d 2 



( e i ' e i) 7^2 + ( e 2 • e 2) irr + (e 3 ■ e 3 



(Aei + Ae 2 + Ae 3 ) . 



OTf OTg ^3 

We can use Clifford to evaluate all these expressions: 
In[l] := « Clifford. m 

The adequate metrics is defined 
In [2]:= $SetSignature = 2 

Here are the differential operators: 
In [3]:= var = {xl,x2,x3}; 

In [4] := GeoGrad[g_] := Sum [InnerProduct [e [k] , e [k] ] *D [g, var [ [k] ] ] *e [k] ,{k,3}] 
In [5]:= GeoDiv[v_]:= Sum [InnerProduct [e [k] , e [k] ] ~2 D [Coef f [v, e [k] ] , var [ [k] ] ] , 

{k,3}] 

In [6] := GeoLap[v_] := Sum [ (InnerProduct [e [k] , e [k] ] ~3) *D [v, {var [ [k] ] ,2}] ,{k,3>] 

The function Coeff [m,b], extracts the coefficient of the blade b in the 
multivector m. The surface and their normal are: 

In [7] := f = xl~2 + x2~2 - x3~2 + R~2; 

In [8]:= norm = GeoGrad[f] / Magnitude [GeoGrad [f] ] ; 

and, finally, the Gaussian curvature is: 

In[9] := kGauss = (InnerProduct [norm, GeoLap [norm] ] +(GeoDiv [norm] ) ~2) /2 

// Simplify- 



Out [9]= 

2 2 2 
xl + x2 - x3 
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Points that lie in the surface fulfills x\ = x\ + x\ + R 2 , therefore: 
In[ll]:= kGauss /. x3 -> Sqrt[xl~2 + x2~2 + R~2] 

-2 

Out [11]= -R 

which is the Gaussian curvature of the hyperboloid. 

Some considerations concerning the isometries (distance-preserving trans- 
formations) of the hyperbolic plane are pertinent before leave this section. 
Like the Euclidean ones, the isometries of the hyperbolic plane can be de- 
scribed in terms of reflections about given axes. All these isometries have 
simple expressions in terms of Clifford algebra. For instance, Eqn. ( 1171) for 
rotations in a given plane remain valid in K 2j i but now is called a Lorentz 
transformation. In general, given the vectors ui, U2, . . . , u& in M p , g , such that 
(ui, u 2 , . . . , Ufc)(ufc, Ufc_i, . . . , ui) = 1, the transformation 

v i — ► (-l) fc (u fc u fc _! • ■ ■ u 1 )v(u 1 u 2 • • - u fe ) = (-l) k UvU, 

is an isometry [l2|. If k is even, the isometry is a rotation, if k is odd, it is 
a reflection [3j. Now, if q = transformations such as the previous one are 
orthogonal. For q = 1 (as the case of the Minkoswky space) they are Lorentz 
t r ansf or mat ions . 

10. Dirac and Pauli algebras 

Now let us add one extra dimension to the space of the previous example 
and consider R3 1. This vector space is fifteen- dimensional with basis 

{1, ei, e 2 , e 3 , e 4 , eie 2 , eie 3 , eie 4 , e 2 e 3 , e 2 e 4 , e 3 e 4 , eie 2 e 3 , eie 3 e 4 , e 2 e 3 e 4 , eie 2 e 3 e 4 }. 

The basis vectors satisfy the relations: 

2_2_2_i 2 _ -i 

e± — e 2 — e 3 — 1, e 4 — —1, 

which is the algebra of the Dirac matrices in relativistic quantum theory. 
Due to this isomorphism, M 31 is often referred as the Dirac algebra. 

The even subalgebra R3 x has a basis {1, eie 2 , eie 3 , eie 4 , e 2 e 3 , e 2 e 4 , e 3 e 4 , e!e 2 e 3 e 4 } 
and is equivalent to the algebra of the Pauli matrices used in quantum me- 
chanics of spin-i particles. This can be seen with the identifications 

o\ = eie 4 , a 2 = e 2 e 4 , cr 3 = e 3 e 4 . 
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By taking i = eie 2 e 3 e4 = criO" 2 o" 3 , we get 

<j\ = 1, i = 1,2,3 
cr l0 - 2 = ia 3 

a 3 ai = ia 2 , 

which are the familiar Pauli matrix relations. is also called the Pauli 
algebra. 

Following the same reasoning, we can prove that the even subalgebra of 
the Pauli algebra is isomorphic to the quaternions. The even subalgebra of 
the quaternions is isomorphic to the complex numbers. The even subalgebra 
of the complex numbers is M.. 



11. Summary 

The basic ideas of the Clifford algebra of a vector space are presented 
and a Mathematica package for calculations within this algebra is developed 
and used in some examples. The relevance of Clifford algebra in physics 
and mathematics lies in the fact that it provides a complete algebraic ex- 
pression of geometric concepts such as directed lines, areas, volumes, etc. 



(For this reason, Clifford algebra is also referred as Geometric algebra [121] ). 
Quantities such as vectors, complex numbers, quaternions, Pauli and Dirac 
matrices, have been normally described by physicists with a potpourri of dis- 
joint mathematical systems. All of them are naturally contained in a Clifford 
algebra. It becomes therefore an efficient mathematical language in a vast 
domain of physics. 
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